Analysis of the rare semileptonic B c — > P{D 1 D s )l + l jvv 
decays within QCD sum rules 



K. Azizi * 

Department of Physics, Middle East Technical University, 06531 Ankara, Turkey 

R. Khosravi ' 

Physics Department , Shiraz University, Shiraz 71454, Iran 



Abstract 

Considering the gluon condensate corrections, the form factors relevant to the 
semileptonic rare B c — ► D,D S (J P = 0~)Z + Z - with I = r, fi,e and B c — > D,D S (J P = 
0~)vis transitions are calculated in the framework of the three point QCD sum rules. 
The heavy quark effective theory limit of the form factors are computed. The branch- 
ing fraction of these decays are also evaluated and compared with the predictions of 
the relativistic constituent quark model. Analyzing of such type transitions could give 
useful information about the strong interactions inside the pseudoscalar D s meson and 
its structure. 
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1 Introduction 



With the chances that in the future a large amount of Be mesons will be produced at LHC 
(with the luminosity values of C = 10 34 cm~ 2 s _1 and y/s = 14TeV, the number of Bf mesons 
is expected to be about 10 s ~ 10 10 per year [1,2]), one might explore the rare Be decays 
to pseudoscalar (D,D S ) and Such types transitions could be useful because of 

the following reasons: 1) Analyzing of such type transitions could give valuable information 
about the nature of the pseudoscalar D s meson and the strong interactions inside it. 2) The 
form factors of these transitions could be used in the study of the polarization asymmetries, 
CP and T violations. 3) These will provide a new framework for more precise calculation 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V tq (q = d, s, b) and leptonic 
decay constants of D sd and B c mesons. 4) These transitions occur at loop level in standard 
model (SM) via the flavour changing neutral current (FCNC) transitions of b — > s, d, which 
are sensitive to the new physics beyond the SM, so these decays are useful to constrain the 
parameters beyond the SM. 5) A possible forth generation, SUSY particles [3] and light 
dark matter [4] might contribute to the loop transitions of b — > s, d. 

The B c , is the only meson containing two heavy quarks with different charge and flavours 
and it is the lowest bound state of b and c quarks, so its decay modes properties are expected 
to be different than flavour neutral mesons. Since the excited levels of be lie below the 
threshold of decay into the pair of heavy B and D mesons, such states decay weakly and 
they have no annihilation decay modes due to the electromagnetic and strong interactions 
(for more about the physics of the B c meson see for example [5]). This paper describes 
the annihilation of the B c into the pseudoscalar (D, D s )l + l~ jvv in the framework of the 
three point QCD sum rules as a non-perturbative approach based on the fundamental QCD 
Lagrangian. This transitions are parameterized in terms of some form factors calculation 
of which plays crucial role in the analyzing of those decay channels. These decays at quark 
level proceed by the loop b — > s, d in the SM with the intermediate u, c and t quarks 
and the main contribution comes from the intermediate top quark. These decay modes 
have also been studied in the relativistic constituent quark model (RCQM) [6]. Some other 
possible channels such as B c — > W'y, B c — > p + 7, B c — > K* + ^, B c — > £?*/ + Z~, B c — > 5*7, 
B c — > D* s d j, B c — ► D* sd l + l~ and B c — »■ Xvv with X be axial vector particle, D s i(2460), 
and vector particles, D*,D* are studied in the light cone or traditional QCD sum rules 
methods in [7-13], respectively. For a set of exclusive nonleptonic and semileptonic decays 
of the B c meson, which have been studied in the relativistic constituent quark model see 
[14]- 

The content of paper is as follows: In section 2, we calculate the sum rules for the related 
form factors considering the gluon correction contributions to the corrolation function. The 
light quark condensate contributions are killed applying the double borel transformations 
with respect to momentum of the initial and final states. The heavy quark effective theory 
(HQET) limit of the form factors are presented in section 3. Section 4 depicts our numerical 
analysis of the form factors and their comparison with the HQET limit of them, results, 
discussions and comparison of our results with the prediction of the RCQM model. 
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Figure 1: loop diagrams for B c — > (D,D s )l + l jvv transitions, bare loop (diagram a) and 
light quark condensates (without any gluon diagram b and with one gluon emission diagrams 
c,d) 

2 QCD Sum rules for transition form factors of the 

B c —> (D, D s )l + l~ I vv 

At quark level, the processes B — > Pl + l~ jvv(P = D,D S ) are described by the loop b — > 
qj. + l~ jvv transitions, (q 1 = d , q 2 = s) in the SM (see Fig.l), and receive contributions from 
photon and Z-penguin and box diagrams for and only Z-penguin and box diagrams 
for vv. These loop transitions occur via the intermediate u, c, t quarks , where dominant 
contribution comes from intermediate top quark. The effective Hamiltonian responsible for 
b — > qj + l~ decays is described in terms of the Wilson coefficients, C^, Cg^ and C w as 



n 



eff 



C e 9 ff qM 1 ' 75)^7/ + C w q^{l - l5 )b£ W i 



2C e 7 ff ^q i ia^(l + l5 )b£ ll ,£ 



(1) 



where Gf is the Fermi constant, a is the fine structure constant at Z mass scale, and Vij 
are elements of the CKM matrix. For vv case, only term containing C±o is considered. The 
amplitudes for for B c — > Pl + l~ jvv decays are obtained by sandwiching of Eq. (1) between 
initial and final meson states: 



C 9 e// < P{p') | qM l - 7 5 )& I B c (p) > I^i 



+ Cio < Pip') I qM l - 7 5 )& I 5 c (p) > i w e 

-2 C 7 e// ^ < P(p') | g- ^"(l + 75)6 I S c (p) > 



(2) 



Next, we calculate the the matrix elements < P(p') | 9^7^(1 — 75)^ | B c (p) > and < 
-P(p') I Qi i^nvQ v {^- + 75)^ I B c (p) > appearing in above equation. The parts of transition 
currents containing 75 don't contribute, so we consider only q. 7^6 and also q i <J tMU q u b parts. 
Considering Lorentz and parity invariances, this matrix elements can be parameterized in 
terms of the form factors as: 



<p{p'MiMBc{p) >-- 

<P(p')\q t ia, u q»b\B c (p) >= 



-(v,U(q 2 )+ q,f~(q 2 )), 

fr(q 2 



[m 



-m P ) 



(3) 
(4) 



m B + Tip 

where f + (q 2 ), f-(q 2 ) and friq 2 ) are the transition form factors, P M = (p + p')^ and q^ = 
{p—p')fi- Here, we should mention that for vv case the form factor fr(q 2 ) does'nt contribute 
since it is related to the photon vertex {a iliV q v ). To calculate the form factors f + (q 2 ), f-(q 2 ) 
and fr(q 2 ), we start with the following correlation function: 



n 



V,T 



where Jp(y) 
of the P and B. 



C759, (g 4 



(5) 



c messons and J„ 



s or d) and J# 

1 - q -,,1' 



= C756 are the interpolating currents 
and Jj = q i iu^ v q u b are transition currents. 
From the general philosophy of the QCD sum rules, we can calculate the above mentioned 
corrolator in two languages: 1) hadron language called the physical or phenomenological 
side, 2) quark gluon language which is the QCD or theoretical side. Equating two sides 
and applying the double Borel transformations with respect to the momentum of the initial 
and final states to suppress the contribution of the higher states and continuum, we get 
sum rules expressions for our form factors. The phenomenological part can be obtained by 
inserting the complete set of intermediate states with the same quantum numbers as the 
currents Jp and Jp c - As a result of this procedure 



(0\Jp\P(p')}(P(p' 



J, 



V.T 



B c {p))(B c {p) 



J 







2 I2\( 2 2\ 

rrip-p )(m Bc -p ) 



(6) 



is obtained. The following matrix elements are defined in terms of the leptonic decay 



constants of the P and B c mesons as: 



(0\J P \P) = -i 



(0\J Bc \B c ) = -i 



fpm 2 p 
m c + m qt 
■ fB c m 2 Bc 



m b + m c 

Using Eqs. (3), (4) and (7) in Eq. (6), we obtain 

fB c m 2 Bc f P m 2 P 



nJ(pW) = 



(m b + m c )(m c + m qi ) (m P - p )(m B -p 



(7) 



f + V, + f.q, +... , (8) 



3 



fB c ml 



f P m 2 P 



(m b + m c )(m c + m qi ) (m 2 P - p )(m 2 B - p 2 ) 



x 



h 



(m Bc + m P ) 



[q 2 V^ - (m| c - m 2 P )qJ\ 



+ ... . 



(9) 



For extracting the expressions for form factors f+{q 2 ) and f-{q 2 ), we choose the coefficients 
of the structures and q^ from UY(p 2 ,p' 2 , q 2 ), respectively and the structure g M from 
Uj^(p 2 ,p' 2 ,q 2 ) is considered for the form factor /t(<? 2 )- Therefore, the correlation functions 
are written in terms of the selected structures as: 



n>V 2 ,g 2 ) = n + ^ + n_g M + 



(10) 



n>W) = Ii T q, + ... . (11) 

On the other side, to calculate the QCD part of correlation function, we evaluate the 
three-point correlator by the help of the operator product expansion (OPE) in the deep 
Euclidean region, where p 2 <C (mj + m c ) 2 and p' 2 <C (m c + m qi ) 2 . For this aim, we write 
each Ilj function in terms of the perturbative and non-perturbative parts as: 



U t (p 2 ,p' 2 , q 2 ) = Ur(pi,pi q 2 ) + Ur per (p\p'\ q 2 ) , (12) 

where % stands for +, — and T and non-perturbative part contains the light quark (< qq >) 
and gluon (< G 2 >) condensates. For the perturbative part, the bare loop diagram (Fig. 1 
a) is considered, however, diagrams b, c, d in Fig. 1 are correspond to the light quark con- 
densates contributing to the correlator. In principle, the light quark condensate diagrams 
give contributions to the correlation function, but applying double Borel transformations 
omits their contributions, hence as first non-perturbative correction, we consider the gluon 
condensate diagrams (see Fig. 2 a, b, c, d, e, f). 

By the help of the double dispersion representation, the bare-loop contribution is written 

as 

lifer _ _ — f f ^ S _P± — ( ' ; Q ) _|_ subtraction terms , (13) 
(2ti) 2 J J (s — p 2 )(s' — p' 2 ) 

where Q 2 = —q 2 . The spectral densities Pi' 6r (s, s', Q 2 ) are calculated by the help of the 
Gutkovsky rule, i.e., the propagators are replaced by Dirac-delta functions 

! -2m5{p 2 - m 2 ) , (14) 



p2 _ m 2 



expressing that all quarks are real. The integration region in Eq. (13) is obtained by 
requiring that the argument of three delta vanish, simultaneously. This condition results 
in the following inequality 

2ss ' + (s + s' + Q 2 )(m 2 -m 2 c -s) + 2s(m 2 c - mgj 

X 1 / 2 (s,s',-Q 2 )X 1 / 2 (mlm 2 c ,s) ' 1 ' 
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where A (a, b, c) = a 2 + b 2 + c 2 — 2ab — 2ac — 2bc. From this inequality, to use in the lower 
and upper limit of the integration over s in subtractions, it is easy to express s in terms of 
s ' i- e - f±( s> ) i n the s — s' plane. 

Straightforward calculations end up in the following results for the spectral densities: 

pl(s,s',q 2 ) = I Nc{A + A' + 

-2m c [(+2 + E 1 + E 2 )m c - (1 + E x + E 2 )m q .\ 
+2m b [(l + E 1 + E 2 )m c - {E 1 + E 2 )m q .} 
+ {E 1 + E 2 )u}, 

pV( s ,s',q 2 ) = I Nc {—A + A' 

-2 m c [(E 2 -E x - l)m q . + (E, - E 2 )m c ] 
-2m b [(l -E x + E 2 )m c +(E 1 - E 2 )m q .] 
+(E 1 - E 2 )u}, 

PT(s,s',q 2 ) = -I N C {A (2m c - m b - m q .) + A' [m b - 2m c + m q .) 
+2[m c (E! - E 2 - 1) + m qi (E 2 - E^} s 
-2 [m 6 (^i - E 2 ) - m c {E x - E 2 + l)]s' 
+ (Ei - E 2 )(m b - 2m c + m q .)u}, 



where 



1 



/o(S ' S '' g) " 4\V 2 ( S , S i,Q 2 Y 
\(s,s',Q 2 ) = s 2 + s' 2 + Q A + 2sQ 2 + 2s Q 2 - 2ss' , 



u 



s + s' + Q 2 , 

2 



A = s + m 2 — m 



a / / i 2 2 

A = s + m, 



c ™,, 



(16) 



and N c = 3 is the color factor. 

Now as first correction to the non-perturbative part of the corrolator, we calculate the 
gluon condensate contributions (see diagrams in Fig. 2). The calculations proceed the 
same as [13] ( see also [9,11,12,18]) and the Fock-Schwinger fixed-point gauge [15-17], 
x^G^ = 0, where is the gluon field is used. In calculations, the following type of 
integrals are encountered: 

d 4 k 1 



I \a, b, c] = [ - — — rr 

J {2tiY [k 2 - m 2 ] a \{p + k) 2 - m 2 ] b \{p' + k) 2 



m 2 
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Figure 2: Gluon condensate contributions to B c — > (D,D s )l + l jvv transitions 



J M [a,6,c] 



d 4 k 



k,. 



(27r) 4 [ k 2 _ m 2^ [(p + fc )2 _ m 2]& (p/ + fc )2 _ m 2 



(17) 



Performing integration over loop momentum and applying double Borel transformations 
with respect to the p 2 and p' 2 , we obtain the Borel transformed form of the integrals as 
follows: 

l 1 \a+b+c 

■(M 2 ) 2 - a - b (M 2 ) 2 - a - c U (a + b + c - 4, 1 - c - b) , 



I (a,b,c) = 



l6ir 2 T(a)T(b)T(c) 



J„(a, b, c) = ^ A(a, 6, c) + f 2 (a, &, c) 



1 



Ji(a, 6, c) - J 2 (a, b, c) 



9m 



(18) 



where 



( _r\a+b+c+l 

h ^ = ' i6^r( a) r( 6 )r( c ) < M ?> 2 '°' t( ^> 3 ""^<'' + » + c - 5, i - « - ») , 

( _l)a+b+c+l 

ls(a ' 6 ' C) = 2 167rW(a)r(&)r(c) (Ml2)3 ' a ' b(M22)2 ' a ' C ^° (a + b + C " 5 ' 1 " C " h) ' 

(19) 

Hat in Eq. (18) denotes the double Borel transformed form of integrals. M\ and Mf are 
the Borel parameters in the s and s' channels, respectively, and the function Uo(a,b) is 
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defined as: 



Wo (a, b) = / dy(y + M 2 + M 2 2 ) a y b exp 
Jo 



-B - B lV 



where 



B 
B l 



1 



.Uf.U.T 

1 



m^Mf + m 2 M 2 4 + M 2 M 2 (m 2 b + m 2 + Q 2 ) 



M{M£ 



[m 2 a + m 2 c )M 2 + M 2 (m 2 b + m 



M 2 M 2 2 ' 



(20) 



After straightforward but lengthy calculations, we get the following results for the gluon 
condensate contributions: 



nf 2 > 



(21) 



where the explicit expressions for Cj are given in appendix-A. 

Next step is to apply the Borel transformations with respect to the p 2 (p 2 — > M 2 ) and 
p/2 _^ on ph enom eno logical as well as the perturbative parts of the correlation 
function, continuum subtraction and equate these two representations of the correlator. 
The following sum rules for the form factors /+ , /_ and fx are derived: 

j + = ( m b + m c ){m c + m qi ) cm 2 n j M?cml/M 2 
+ fB c m 2 B J P m 2 P 



x 



47T 2 I J(m c +rn q )' 2 



ds' / 

Jf-W) 



min(s , /+(«')) 



dspl(s, s', Q 2 )e- s ' M h- s 'l M l - %M\Ml ( —G 

\ 7T 



c 4 



(22) 



j_ = ( m b + m c )(m c + "O e m| e /M? e m»/M3 

fB c m 2 B J P m 2 P 



X 



ds' 



47T 2 I J(m c +m q ) 2 J /_(«') 



min(s J + (s')) 



dsp v (s, s', Q 2 )e- S / M h- S '/ M l - %M\Ml l^G 2 ) — 

\ 7T / 6 



(23) 



y T = (m b + mj K + m 9i ) ^ /M 2 ^ /M 2 
T fB c m 2 B J P m 2 P (m Bc - m P ) 

47T 2 \ J(m c +m q )l Jf.(s') ' t> / 1 2 \ 7T / 6 J 

(24) 

where So and s' are the continuum thresholds and s = f±(s') in the lower and upper limit 
of the integral over s are obtained from inequality ( 15). The min(so, f+(s')) means that for 
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each value of the q 2 , the smaller one between sq and /+ is selected. In above equations, in 
order to subtract the contributions of the higher states and the continuum the quark-hadron 
duality assumption is also used 



p hi 9 herstates^ ^ = p OPE^ _ _ 



(25) 



At the end of this section, we would like to present the differential decay width of 
B c — > Pl + l~ jvv decays. Using the parametrization of these transitions in terms of form 
factors and amplitude in Eq.(2), we get 



dT 



B 



± 



P ± vu) 



v tq y; b \ 2 ( t>T{i,r p ,s)m i Bc \c 



dQ 2X c ' 2 8 7T 5 

where, 0p(l,r p , s) is the usual triangle function 



10 



2\|2 



(26) 



;i, r p , s) = 1 + r 2 p + s 2 - 2r p - 2s - 2r p s, 



with r 



m 



Q 2 



P 2 ' 

m B c 



2 ' 
m B c 



and 
dT 



dQ 



2 (Bt^p±in-) = 



G 2 F \V tq V t * b \ 2 ml a 2 



3 • 2 9 7T 5 



v(f> p (l,r P ,s) 



s 



'l,r P ,s)a 1 + 12tp 1 

(27) 



where t = mf /m 2 3c and the expressions of a± and (5\ and v are given as: 



1 + 



4m? 



2m b Cf f T (Q 2 



1 + 



+ \c w f+(Q 



2\|2 



A = |C : 



10 | 



1 



l2\|2 



1 + r P - - + 1 - r P )Re(f + (Q 2 )f_(Q 2 )) + -s\f_(Q 2 ) 



where m b = mb/rriB c - 



3 HQET limit of the form factors 

In this section, we present the infinite heavy quark mass limit of the form factors for 
B c — > (D, D s )l + l~ jvv transitions. To this aim, we use the following parametrization (see 
also [19-23]): 

m r + m 2 n — q 2 

y = vv' = (28) 

2m B m Dqi 

where v and v 1 are the four- velocities of the initial and final meson states, respectively and 
y = 1 are so called zero recoil limit. Now, to obtain the y dependent expressions of the 
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form factors we define m& — > oo, m c = where z is given by \fz = y + — 1 and we 
also set the mass of light quarks to zero. In this limit the new Borel parameters 7\ and T 2 
take the form T\ = Mf/2m b and T 2 = M 2 2 /2m c . 

The new continuum thresholds z/ , and z/q are defined as: 

Vo = ftZlS, K = fo^I, (29) 
m b m c 



and the new integration variables become 



v = 1/ = (30) 

The leptonic decay constants are rescaled: 

he = V™bfB c , /d,. = V™~cfD q .- (31) 

The correspond expressions for I (a,b,c), Ii(a,b,c) and I 2 (a,b,c) in this limit are defined 
as: 

( 1 \a+ft+c 

C)HQET = leJrwrwrM ^' 2 ""''^' 2 '"^" "'" + » + «- 4,i-c-»), 

/_i \a+b+c+l 

h[a ' "• C) " QET = ' i6^r(a)r(6)r ( e) (r ' )2 ' a ' t(:r2)3 '°'^" OET(a + " + c - 5, 1 - c - 6) , 

/i \a+6+c+l 

^ * C) "° £r ^ + » + c - 5, 1 - c - ») , 

(32) 

where 7\ and T 2 are the Borel parameters in the s and s' channel, respectively, and the 
function Uq QE1 \m,n) is defined as 

Uo QET {m, n) = f°°(x + Ti + T 2 ) m :r n exp{i + S + C}dx, (33) 
Jo 

with 

m h 2 T 2 2 + T x T 2 {m h 2 + y) 



A = 
B = 



T 1 T 2 x 
T 2 (m b 2 + + 
T X T 2 



C = (34) 

In order to the calculations be easy, the following redefinitions for the form factors are 
applied 

fi = fi{m Bc + m Dq .} (35) 
9 



After some calculations, we obtain the y-dependent expressions of the form factors as 
follows: 



7HQET 
J+ 



(y) 



1 



a if 3(1 + V /I) 5 



^ 2 f B J p 



f -l + y/Z)F(y,z) 



[-1 + (1 + 3y)^ - (2 + y + Ay 2 )z + (3 + 2y(2 + y))z 3/2 - Ayz 2 + 2z 5/2 ] 



J du J due 2T ie 2T ^Q{2yvv 



(36) 



1 {V> ZWfBjp \(-l + y/z)F(y,z) 



[-1 + (1 + 3y)v^ - (2 + 5y)^ + (7 + 2y(y - 2))z 3/2 - 4(y - l)z 2 + 2z 5/2 \ 

/ du due 2T ie 2T 20(2yuu - u 2 - u' 2 ) 
Jo Jo 



= L_e^e*/ -3(1 + v^) 



where 



32n*f B J P [(-l + y/Z)F(y,z) 

[3 - (1 + %)v^ + (4 + y(3 + 8y))z - ((2 + y)(l + Ay))z 3 ' 2 + (3 + 2y(y + 2))z 2 

-Ayz 5/2 + z 3 ] [ ° du [ ° du'e~^e~^9(2yuu - u 2 - u' 2 ) 
Jo Jo 

+ hm (i -J—, (1 + v^) 2 / ^G 2 ) C« QET ) ) , (38) 
F(y, ^) = ^ 3/4 [l + z + y 2 z + 2 2 - 2^(1 + ^)] 3/2 . (39) 



In the heavy quark limit expressions of the form factors, A = va Bq — and A = mo* — ni c 
and the explicit expressions of the coefficients C^ ET are given in the appendix-B. 



4 Numerical analysis 

This section encompasses our numerical analysis of the form factors /+ , /_ and f T and 
their HQET limit, branching fractions, comparison of our results with the prediction of the 
RCQM and discussion. The sum rules expressions of the form factors depict that the main 
input parameters entering the expressions are gluon condensate, Wilson coefficients C? , 
Cg and Cio , elements of the CKM matrix V t b, V ts and V t d, leptonic decay constants; 
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fB C i Jd and fo s , Borel parameters M 2 and Mf, as well as the continuum thresholds so 
and s' . In further numerical analysis, we choose the values of the Gluon condensate, 
leptonic decay constants, CKM matrix elements, Wilson coefficients, quark and meson 
masses as: < ^G 2 >= 0.012 GeV 4 [24], C e 7 ff = -0.313, C 9 e// = 4.344, C w = -4.669 
[25,26], | V tb |= 0.77tg;& | V u \= (40.6 ± 2.7) x 10~ 3 | V td \= (7.4 ± 0.8) x 10~ 3 [27], 
f Ds = 274±13±7 MeV [28], f D = 222.6 ± 16.7±1;| MeV, [29], f Bc = 350±25 MeV [30-32], 
m c = 1.25 ± 0.09 GeV, m s = 95 ± 25 MeV , m fe = (4.7 ± 0.07) GeV, m d = (3 - 7) MeV, 
m Ds = 1.968 GeF, m D = 1.869 GeF, m Bc = 6.258 GeV [33], A = 0.62GeV [34] and 
A = 0.86GeV[35]. 

The expressions for the form factors contain also four auxiliary parameters: Borel mass 
squares Mf and Mf and continuum threshold sq and s' . These are not physical quantities, 
so the physical quantities, form factors, should be independent of them. The parameters 
so and s' , which are the continuum thresholds of B c and P mesons, respectively, are deter- 
mined from the conditions that guarantees the sum rules to have the best stability in the 
allowed M\ and Mf region. The values of continuum thresholds calculated from the two- 
point QCD sum rules are taken to be s = (45 - 50) GeV 2 and s' = (6 - 8) GeV 2 [7, 24, 36] . 
The working regions for M\ and M| are determined by requiring that not only contribu- 
tions of the higher states and continuum are effectively suppressed, but the gluon condensate 
contributions are small, which guarantees that the contributions of higher dimensional op- 
erators are small. Both conditions are satisfied in the regions 10 GeV 2 < Mf < 25 GeV 2 
and 4 GeV 2 < Mf < 10 GeV 2 . 

The dependence of the form factors /+ , /_ and fx on M 2 and Mf for B c — > D s l + l~ jvv 
are shown in Figs. 3, 4 and 5, respectively. The figures 6, 7, and 8 also depict the dependence 
of the the form factors on Borel mass parameters for for B c — > Dl + l~ jvv. This figures show 
a good stability of the form factors with respect to the Borel mass parameters in the working 
regions. Our numerical analysis shows that the contribution of the non-perturbative part 
(the gluon condensate diagrams ) is about 8°/o of the total and the main contribution comes 
from the perturbative part of the form factors. 

The values of the form factors at q 2 = are shown in Table 1: The sum rules for the 







B C ^D S 


U{in-/w) 


0.22 ± 0.045 


0.16 ±0.032 


f-{l + l-/uu) 


-0.29 ±0.056 


-0.18 ±0.038 


fr{l + l-) 


-0.27 ±0.054 


-0.19 ±0.040 



Table 1: The values of the form factors at q 2 = 

form factors are truncated at about 2 GeV 2 below the perturbative cut, so to extend our 
results to the full physical region, we look for parametrization of the form factors in such 
a way that in the region < q 2 < 19.26 (18.41) GeV 2 for D(D S ), this parametrization co- 
incides with the sum rules prediction. Our numerical calculations shows that the sufficient 
parametrization of the form factors with respect to q 2 is as follows: 

fM 2 ) = nr^x^ ( 4 °) 

1 ± aq + pq l 
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where q = q 2 /m 2 Bc . The values of the parameters fi(0), a and (3 are given in the Tables 2, 3. 





f(0) 


a 


P 


U{in-/uu) 


0.22 


-1.10 


-2.48 


f-{in-/vD) 


-0.29 


-0.63 


-4.06 


f T {i + n 


-0.27 


-0.72 


-3.24 



Table 2: Parameters appearing in the form factors of the B c — > Dl + l jvv decay for M\ = 
15 GeV 2 , Ml = 8 GeV 2 . 





f(0) 


a 


P 


U{in-/uu) 


0.16 


-1.55 


-2.80 


f-{in-/vu) 


-0.18 


-0.77 


-6.71 


f T (i + n 


-0.19 


-1.43 


-3.06 



Table 3: Parameters appearing in the form factors of the B c — > D s l + l jvv decay for 
M\ = 15 GeV 2 , Ml = 8 GeV 2 . 

The errors are estimated by the variation of the Borel parameters M\ and Mf , the variation 
of the continuum thresholds sq and s' Q , the variation of b and c quark masses and leptonic 
decay constants Jb c and fo,(D s )- The main uncertainty comes from the thresholds and the 
decay constants, which is about ~ 18% of the central value, while the other uncertainties 
are small, constituting a few percent. 

Now, we compare the extrapolation values for the form factors and their HQET values 
obtained from Eqs. (36-38) in Tables 4 and 5 for B c — > Dl + l~/vv and B c — > D s l + l~ /vv, 
respectively. 



y 


i 


1.1 


1.2 


1.3 


1.4 


1.5 


1.6 


1.7 


1.8 


g 2 


19.26 


16.93 


14.59 


12.25 


9.91 


7.57 


5.23 


2.89 


0.55 


UW Z ) 


2.19 


1.36 


0.88 


0.56 


0.36 


0.29 


0.27 


0.24 


0.23 


f-W 2 ) 


-3.01 


-1.93 


-1.20 


-0.75 


-0.52 


-0.39 


-0.33 


-0.32 


-0.31 




-2.52 


-1.53 


-1.12 


-0.70 


-0.49 


-0.37 


-0.31 


-0.29 


-0.28 


fi Qm \y) 


? 


1.35 


0.50 


0.29 


0.20 


0.15 


0.12 


0.10 


0.08 




? 


-1.90 


-0.75 


-0.44 


-0.30 


-0.22 


-0.18 


-0.15 


-0.12 


f? QET (y) 


? 


-1.51 


-0.58 


-0.33 


-0.23 


-0.17 


-0.14 


-0.11 


-0.10 



Table 4: The comparison of the extrapolation values for the form factors and their 
HQET limit for B c -> Dl + l~ at M 2 = 15 GeV 2 , Ml = 8 GeV 2 and corresponding 
Ti = 1.6 GeV,T 2 = 3.2 GeV. 

At the y — 1 called the zero recoil limit, the HQET limit of the form factors are not 
finite and at this value, we can determine only the ratio of the form factors. For other 
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y 


1 

1 


i i 
i . i 


1 9 


1 3 

l.O 


1 zL 


1 ^ 

l.O 


l.U 


1 7 


a' 2 


1 % zLI 


1 ^ QzL 


10.40 




o.oo 




o.uo 


1 1 fi 


f , (a 2 ) 


2 1 7 


1.12 


79 


53 


31 


0.22 


18 


17 


J — \H ) 


-2.50 


-1.53 


-0.79 


-0.43 


-0.29 


-0.24 


-0.23 


-0.22 


fM 2 ) 


-2.25 


- 1.23 


-0.70 


-0.37 


-0.27 


-0.23 


-0.21 


-0.20 


fH Qm Xv) 


? 


1.08 


0.41 


0.24 


0.16 


0.12 


0.10 


0.08 


f HQE1 \y) 


? 


-1.52 


-0.60 


-0.35 


-0.24 


-0.18 


-0.14 


-0.12 


fr K1 (y) 


? 


-1.22 


-0.46 


-0.27 


-0.18 


-0.14 


-0.11 


-0.10 



Table 5: The comparison of the extrapolation values for the form factors and their 
HQET limit for B c -> DJ+l' at M 2 = 15 GeV 2 , M 2 = 8 GeV 2 and corresponding 
T x = 1.6 GeV,T 2 = 3.2 GeV. 

values of y and corresponding q 2 , the behavior of the form factors and their HQET values 
are the same, i.e., when y increases (q 2 decreases) both the form factors and their HQET 
values decrease. Moreover, at high q 2 values, the form factors and their HQET values are 
close to each other while at low q 2 , the form factor values are about 2-3 times greater than 
that of their HQET limit. 

At the end of this section we would like to present the values of the branching ratios. 
Integrating Eqs. (26) and (27) over q 2 in the whole physical region and using the total mean 
life time r ~ 0.46 ps of B c meson [37], the branching ratio of the B c — > P(D, D s )l + l~ jvv 
decays are obtained as Table 6. This Table also includes a comparison of our results 
with the prediction of the RCQM. This Table presents a good agreement between two 
models especially when the errors are taken into account. Any experimentally measurements 
on the branching fractions of these decays and those comparisons with the results of the 
phenomenological models like QCD sum rules could give valuable information about the 
nature of the D s meson and strong interactions inside it. 

In summary, we investigated the B c — > P(D, D s )l + l~ jvv channels and computed the 
relevant form factors and their HQET limits considering the gluon condensate corrections. 
We also evaluated the total decay width and the branching fractions of those decays and 
compared our results with the predictions of the RCQM. Detection of these channels and 
their comparison with the phenomenological models like QCD sum rules could give useful 
information about the structure of the D s meson. 
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decay 


our results 


RCQM[6] 


JD C ? U IS IS 


^ 48 -1- 71 "\ v 1 n~ 8 


O.^O A IU 


R — > D in/ 

J->C LJ gIS IS 




U.I A 1 U 


r — i- r>p+p- 


h ^4 4- n 9^ v i n~ 8 






/i A 7 _|_ n oo\ v i n-7 
i ± i zn u.Oi i a iu 




r — r>ii + n~ 


l u.oi zn U.UU 1 A _LU 


n 44 v 1 n~ 8 


B c —> D s fi + fi~ 


(0.61 ±0.15) x 10 _Y 


0.97 x 10" Y 


B c -> Dr+r- 


(0.13 ±0.03) x 10" s 


0.11 x 10~ s 




(0.23 ±0.05) x 10"'' 


0.22 x 10- y 



Table 6: Values for the branching fractions of the B c — > P(D, D s )l + l jvv decays and their 
comparison with the predictions of the RCQM [6] 
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Appendix-A 

In this appendix, the explicit expressions of the coefficients of the gluon condensate 
entering the sum rules of the form factors /+, /_ and fx are given. 



C+ = -5/ 1 (3,2,2)m c 6 -5/ 2 (3,2,2)m c 6 -5/ (3,2,2)m c 6 
±5 7 2 (3, 2, 2)m c 5 m b + 5 h(3, 2, 2)m c 5 m b ± 5 7 (3, 2, 2)m c 5 m h 
±5 7 2 (3, 2, 2)m c 4 m b 2 + 5 h(3, 2, 2)m c A m b 2 - 5 ii(3, 2, 2)m c 3 m 6 3 
-5 7 (3, 2, 2)m c 3 m 5 3 - 5 J 2 (3, 2, 2)m c 3 m 6 3 - 5 7 2 (3, 2, l)m c 4 
±15 ij ' 1] (3, 2, 2)m c 4 - 15 h(2, 2, 2)m c 4 ± 15 7 2 ' 1] (3, 2, 2)m c 4 
-5 J 2 (3, 1, 2)m c 4 - 5 h(3, 1, 2)m c 4 - 15 J 2 (4, 1, l)m c 4 
-5 7^3, 2, l)m c 4 - 15 7 (4, 1, l)m c 4 - 15 7 (2, 2, 2)m c 4 
±15 7i°' 1] (3,2,2)m c 4 - 15 7^4, 1, l)m c 4 - 15 7 2 (2, 2, 2)m c 4 
±5 7i(3, l,2)m c 3 m 6 - 10 (3,2,2) m c 3 m b + 10 7i(2, 2, 2)m c 3 m 6 
±10 7 (3, 2, l)m c 3 m 6 - 10 7f ) ' 1] (3, 2, 2)m c 3 m 6 + 5 7 (3, 1, 2)m c 3 m 6 
-10 7 2 (2, 3, l)m 3 m b + 15 7 (4, 1, l)m c 3 m 6 ± 5 7i(3, 2, l)m c 3 m£, 
±15 7i(4, 1, l)m 3 m b - 10 7^2,3, l)m 3 m b + 15 7 2 (4, 1, l)m 3 m b 
±10 7 2 (2, 2, 2)m c 3 m 6 - 10 7 (2, 3, l)m 3 m b - 10 7{°' 1] (3, 2, 2)m c 3 m fe 
±5 7 2 (3, 1, 2)m 3 m b + 10 7 (2, 2, 2)m 3 m b + 5 7 2 (3, 2, l)m 3 m b 
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m c z m b 2 
m b 



m c 2 
m r 2 



+5 4°' 1] (3, 2, 2)m c 2 m 6 2 + 30 J„(l, 4, l)m c W - 5 J (2, 2, 2)m c 2 
-10/i(3,2, l)m c 2 m 6 2 + 307 2 (1,4, l)m c 2 m 6 2 - 107 2 (3,2, l)m c 2 
+30 7i(l, 4, l)m c 2 m 6 2 - 5 4°' 1] (3, 2, 2)m c m b z - 30 7 (1, 4, l)m c m b 3 
+ 10 7 2 (2, 3, l)m c m 6 3 + 10 7^2, 3, l)m c m 6 3 - 5 7{°' 1] (3, 2, 2)m c mt, 3 
-5 7 2 °' 1] (3, 2, 2)m c m b 3 + 10 7 2 (3, 2, l)m c m 6 3 - 30 7^1, 4, l)m c m b z 
+ 10 7i(3, 2, l)m c m b 3 - 30 7 2 (1, 4, l)m c m 6 3 + 15 7 (1, 4, l)m 6 4 
-5 7 (3, 2, l)m b A + 15 4°' 1] (3, 2, l)m c 2 - 5 7 (3, 1, l)m c 2 
-5 7 (2, 2, l)m c 2 - 5 7 (2, 1, 2)m c 2 + 30 4°' 1] (2, 2, 2)m c 2 
+ 15 7^(4, 1, l)m c 2 + 15 7 2 °' 1] (3, 2, l)m c 2 + 20 7i°' 11 (3, 2, 1). 
-157f 2l )(3,2,2)m c 2 + 20 7i ' 11 (3,l,2)m c 2 + 154 ' 11 (4,l,l) b 
+ 15 7{°' 1] (4, 1, l)m c 2 - 15 7{°' 2] (3, 2, 2)m c 2 + 15 7 2 °' 1] (3, 1, 2)m c 2 
-10 7 1 (l,2,2)m c 2 - 10 7 (l,2,2)m c 2 + 30 7i 0,11 (2,2,2)m c 2 
+30 (7 2 °' 1] (2,2,2)m c 2 - 10 7 2 (1, 2 , 2)m c 2 + 15 7{°' 1] (3, 1, 2)m c 2 
-15 7i 0,2] (3, 2, 2)m c 2 - 25 7 (2, 2, l)m c m b - 5 7 2 (2, 1, 2)m c m 6 
-10 7[°' 1] (3,2,1) m c m 6 - 20 7i(2, 2, l)m c m 6 - 40 7 (1, 3, l)m c m 6 
-20 7 2 (2,2, l)m c m 6 - 10 7{ 0,1] (2, 2, 2)m c m b - 10 7i(l, 2, 2)m c m 6 
-10 /| 0,1] (2, 2, 2)m c m 5 - 40 7^1, 3, l)m c m 6 - 5 7^2, 1, 2)m c m 6 
-5 7 2 °' 1] (3,2, l)m c m 6 - 10 4° ,1] (2, 2, 2)m c m 6 - 15 7{ 0,1] (3, 1, 2)m c m b 
-5 7 (2, 1 , 2) m c m b + 5 7{ 0,2] (3, 2, 2) m c m 5 - 15 7i° ,1] (3, 1 , 2) m c m b 
+10 7i°' 1] (2,3, l)m c m b - 10 7 2 (1, 2, 2)m c m 6 - 10 7 (1, 2, 2)m c m 6 
-15 7 2 0,1] (3, 1, 2) m c m 6 + 5 7 2 °' 2] (3, 2, 2) m c m b - 5 7} 0,11 (3, 2, 1) m c m b 
+5 7i 0,2] (3, 2, 2)m c m 6 - 40 7 2 (1, 3, l)m c m 6 + 10 7{°' 1] (2, 3, l)m c m 6 
+ 10 7 2 °' 1] (2,3,l)m c m b + 10 7{°' 1] (2, 2, 2)m 6 2 - 30 7 2 °' 1] (1, 4, l)m 6 2 
+ 15 7 (1, 3, l)m 6 2 + 5 7 (2 , 2, l)m 6 2 + 5 7 2 (2 , 2, l)m 6 2 
-5 7 2 °' 2] (3, 2, 2)m 6 2 - 5 7 (3, 1, l)m 6 2 - 5 7 (1, 2, 2)m 6 2 
+ 10 7{ ' 1] (3,2,l)m 6 2 -5 7l ' 2l (3,2,2)m & 2 - 30 7i ' 11 (l,4,l)m b 2 
+ 10 7 2 °' 1] (3, 2, l)m 6 2 + 10 7 2 °' 1] (2, 2, 2)m 6 2 - 30 4°' 1] (1, 4, l)m 6 2 
-5 7 (2, 1, 2)m b 2 + 15 7i°' 1] (2, 2, 2)m b 2 + 5 7i(2, 2, l)m 6 2 
-107{°' 2] (3, 1, 2) + 10 7 2 (1, 2, 1) + 5 7 (2, 1, 1) + 10 7' ' 11 (2, 1, 2) 
+ 10 7x(l, 1, 2) + 10 7 2 (1, 1, 2) + 5 7 (1, 2, 1) + 10 4°' 1] (2, 2, 1) 
-10 7 2 0,2] (3,2, 1) - 15 7 2 0,2] (2,2,2) - 10 7{ 0,2] (3, 2, 1) + 5 7 (1, 1, 2) 
+5 4°' 1] (2, 2, 1) + 5 7 2 °' 1] (2, 1, 2) + 5 7{°' 1] (2, 1, 2) + 5 7 2 °' 1] (2, 2, 1) 
-15 4°' 2] (2, 2, 2) - 15 /{°' 2] (2, 2, 2) + 10 h(l, 2, 1) + 10 h(2, 1, 1) 
+ 10 7 2 (2, 1, 1) + 10 4°' 1] (3, 1, 1) - 10 7 2 °' 2] (3, 1, 2) + 10 7<f ' 1] (1, 2, 2) 
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+10 7} U ' 1J (1,2,2) + 10 7 2 U ' 1J (1,2,2) , 

C_ = 5 7 2 (3, 2, 2)m c 6 - 5 h(3, 2, 2)m c 6 - 5 7 (3, 2, 2)m c 5 m b 
+5 7i(3, 2, 2)m c 5 mf, - 5 7 2 (3, 2, 2)m c 5 m 6 + 5 7i(3, 2, 2)m c 4 m 6 2 
-5 7 2 (3, 2, 2)m c 4 m b 2 - 5 7^3, 2, 2)m c 3 m 6 3 + 5 7 (3, 2, 2)m c 3 m 6 3 
+5 7 2 (3, 2, 2)m c 3 m 5 3 - 15 7 2 °' 1] (3, 2, 2)m c 4 - 5 7i(3, 1, 2)m c 4 
+15 7 2 (2, 2, 2)m c 4 - 15 7^4, 1, l)m c 4 + 5 7 2 (3, 2, l)m c 4 
+5 7 2 (3, 1, 2)m c 4 + 15 7j°' 1] (3, 2, 2)m c 4 - 5 7^3, 2, l)m c 4 
-15 7i(2, 2, 2)m c 4 + 15 7 2 (4, 1, l)m c 4 + 5 7^3, 1, 2)m 3 m b 
+ 10 7 2 0,11 (3, 2, 2)m 3 m b + 10 /i°' 1] (3, 2, 2)m c 3 m 6 - 5 7 2 (3, 1, 2)m 3 m b 
-10 7i(2, 3, l)m 3 m b - 15 7 (4, 1, l)m c 3 m b + 10 7 2 (2, 3, l)m 3 m b 
+15 7^4, l,l)m c 3 m 6 - 15 7 2 (4, 1, l)m 3 m b - 10 7 (2, 3, l)m 3 m b 
+5 7 (3, 2, l)m 3 m b - 5 7 (3, 1, 2)m 3 m b - 10 7 2 (2, 2, 2)m c 3 m 6 
-10 7 (2, 2, 2)m 3 m b - 10 7{ 0,1] (3, 2, 2)m 3 m b + 5 7^3, 2, l)m 3 m b 
-5 7 2 (3, 2, l)m c 3 m 6 + 10 7i(2, 2, 2)m c 3 m 6 - 5 4°' 1] (3, 2, 2)m c 2 m b 2 
+5 7 (2, 2, 2)m c 2 m 6 2 + 30 7^1, 4, l)m c 2 m 6 2 + 10 7 2 (3, 2, l)m c 2 m 6 2 
-10 7i(3, 2, l)m c 2 m b 2 - 30 7 2 (1, 4, l)m c 2 m 6 2 + 5 I [ q' 1] (3, 2, 2)m c m b 3 
+30 7 (1, 4, l)m c m 6 3 - 30 7i(l, 4, l)m c m 6 3 + 10 7i(2, 3, l)m c m b 3 
+10 7i(3, 2, l)m c m b 3 - 5 7{°' 1] (3, 2, 2)m c m b 3 + 30 7 2 (1, 4, l)m c m 6 3 
-10 7 2 (2, 3, l)m c m 6 3 - 10 7 2 (3, 2, l)m c m 5 3 + 5 7 2 °' 1] (3, 2, 2)m c m 6 3 
+5 7 (3, 2, l)m & 4 - 15 7 (1, 4, l)m 6 4 + 15 4° ,1] (4, 1, l)m c 2 
-10 7 (2, 2, l)m c 2 + 10 4°' 1] (3, 2, l)m c 2 + 10 7 (2, 1, 2)m c 2 
+30 7{ 0,11 (2,2,2)m c 2 + 10 7 2 (1, 2 , 2)m c 2 - 10 7i(l, 2, 2)m c 2 
-30 7 2 °' 1] (2, 2, 2)m c 2 - 15 7 2 °' 1] (3, 1, 2)m c 2 - 15 7 2 °' 1] (4, 1, l)m c 2 
+15 7{°' 1] (3, 1, 2)m c 2 - 15 7 2 °' 1] (3, 2, l)m c 2 + 15 7{°' 1] (3, 2, l)m c 2 
-10 4°' 1] (3, 1, 2)m c 2 - 15 7}° ' 2] (3, 2, 2)m c 2 + 15 7 2 °' 2] (3, 2, 2)m c 2 
+40 7 2 (1, 3, l)m c m b - 5 4°' 2] (3, 2, 2)m c m 6 - 5 7{°' 1] (3, 2, l)m c m 6 
+ 15 7i 0,11 (3,l,2)m c m 6 - 10 7^1, 2, 2)m c m 6 + 15 7 2 °' 1] (3, 1, 2)m c m b 
-15 7}°' 1] (3, l,2)m c m b - 10 7}°' 1] (2, 2, 2)m c m b + 5 7 2 0,1] (3, 2, l)m c m 6 
+20 7 (1, 3, l)m c m 6 + 5 7j°' 2] (3, 2, 2)m c m 6 + 5 7 2 (2, 1, 2)m c m 6 
+20 7 (2, 2, l)m c m 6 - 20 7^2, 2, l)m c m 6 + 10 7 2 (1, 2 , 2)m c m b 
+10 7i 0,1] (2, 2, 2)m c m b - 10 7 2 °' 1] (2, 3, l)m c m 6 + 15 4° ,1] (3, 2, l)m c m 6 
-5 7i(2,l, 2) m c m b + 10 7 2 0,1] (2, 2, 2)m c m b + 10 7{ 0,1] (2, 3, l)m c m 6 
+107 (1, 2, 2)m c m 6 + 5 7 (2, 1, 2)m c m 6 - 5 7 2 °' 2] (3,2,2) 
+20 7 2 (2, 2, l)m c m 6 - 40 7^1, 3, l)m c m 6 + 10 4° ,1] (2, 3, l)m c m 6 
-15 7 (1, 3, l)m 6 2 - 30 4, l)m fe 2 + 5 7 (2, 1, 2)m 6 2 
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-10 7 2 °' 1J (3, 2, l)m b 2 + 10 4°' 1| (3, 2, l)m b 2 - 5 J 2 (2, 2, l)m 6 2 
+5 Ji(2, 2, l)m fe 2 - 5 /J 0,21 (3, 2, 2)m 6 2 + 5 J (3, 1, l)m 6 2 
+5 4°' 2] (3, 2, 2)m 6 2 + 30 4°' 1] (1, 4, l)m 6 2 - 5 if ' 1] (2, 2, 2)m 6 2 
+10 4°' 1] (2, 2, 2)m 6 2 - 10 J„(2, 2, l)m 6 2 - 10 4°' 1] (2, 2, 2)m 6 2 
-10/ 2 (2, 1, 1) + 107^1, 1,2) - 10/ 2 (1, 1,2) + 10/f' 1] (l,2,2) 
-10 7 2 (1, 2, 1) + 5 4°' 1] (2, 2, 1) + 10 4°' 2] (3, 2, 1) - 15 /}°' 2] (2, 2, 2) 
+15 4 ' 21 (2, 2, 2) + 5 4°' 1] (2, 1, 2) - 10 /J 0,21 (3, 2, 1) - 10 /} ' 21 (3, 1, 2) 
+10 4°' 2] (3, 1, 2) - 5 7<f' 1] (2, 1, 2) - 10 4°' 1] (1, 2, 2) +5/{°' 1] (2, 2, 1) 
-5 4°' 1] (2, 1, 2) - 5 4°' 1] (2, 2, 1) + 10 h(l, 2, 1) + 10 7 1 (2, 1, 1) . 



where 



C T = -5 7 (3, 2, 2)m c W + 5 J (3, 2, 2)m c W - 5 J (3, 1, 2)m c 4 m fe 
+5 J (3, 2, l)m c 4 m 6 - 5 J (3, 1, 2)m c 3 m b 2 + 5 J (3, 2, l)m c 3 m b 2 
+ 10 4°' 1] (3, 2, 2)m c 2 m b 3 - 15 7 (3, 2, l)m c 2 m b 3 - 15 7 (4, 1, l)m c 2 m 6 3 
-10 7 (2, 3, l)m c 2 m 6 3 + 5 7 (3, 1, 2)m c 2 m 6 3 - 10 7 (2, 2, 2)m c 2 
+5 4°' 1] (3, 2, 2)m b 5 - 5 7 (2, 2, 2)m b 5 - 10 7 (2, 3, l)m b 5 
-107 (3,2, l)m 6 5 + 307 (1,4, l)m b 5 + 10 7 (3, 1, l)m c 3 
-10 7 (2, 1, 2)m c 2 m 6 - 10 7 (3, 1, l)m c 2 m fe - 10 4°' 1] (3, 2, l)m c 2 m b 
+ 10 4°' 1] (3, 1, 2)m c 2 m b - 5 4°' 1] (3, 2, l)m c m b 2 + 5 7 (2, 2, l)m c m 6 2 
-5 7 (2, 1, 2)m c m b 2 + 5 4° ' 1] (3, 1, 2)m c m fc 2 - 10 7 (1, 2, 2)m b 3 
-5 7 (2, 2, l)m 6 3 - 5 4°' 2] (3, 2, 2)m fe 3 + 50 7 (1, 3, l)m 5 3 
+ 15 4° ,1] (3, 1, 2)m b 3 + 5 4°' 1] (3, 2, l)m b 3 - 10 7 (3, 1, l)m fe 3 
-10 7 (2, 1, 2)m fe 3 + 10 7' ' 11 (2, 3, l)m fc 3 + 10 7' ' 11 (2, 2, 2)m 6 3 
— 104°' 1] (3 ! 1, l)™c + 10 7 (2, 1, l)m c + 30 7 (1,2, l)m b 
-10 7 (1, 1, 2)m 6 + 5 7 (2, 1, l)m b - 5 4°' 2] (3, 1, 2)m fe 
+5 7i 0,2l (3,2, l)m b + 10 7' ' 11 (2, l,2)m 6 



7^(a,6 )C ) = (M 2 ) l (M 2 ) J 



d (Mi 2 )' d (M|) J 

Appendix— B 



^[(M 2 ) l (M 2 ) J 7 n (a,fe,c) 



In this appendix, the explicit expressions of the coefficients of the gluon condensate 
entering the HQET limit of the form factors f+ T , f^ ET an d T are given. Note 
that only in this appendix, by Ii(a,b, c) we mean Ii(a, b, c) H ® ET wh.ich are defined in Eq. 
(32). 



gHQET 



4°' 1] (3,2,2)m 6 5 , 7o'' iJ (3,2,2)m fe 4 oo 7 2 (2, 1, l)m b 4 



?[o,i], 



+ 



-32 
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_ 4 J (2,2,l)m fc 4 _ g I (2,l,l)m b 3 + ^ 4°' 2] ( 3, 1, 2)m 6 3 



_ g y'lM.lK + 4 I (2,l,2)m b 4 + 12 J 2 u ' iJ (3, 1, 2)m fe 4 



?[o,i], 



2 v- 2 v- 2 

i6 / (l,l,2)m b 4 _ 12 4°' 1] ( 2, 2, 2)m fe 4 Q / 2 °' 1] (3, 2, l)m 6 4 



+g 7 (l,2,2)m„ 5 + g I 2 (2,l,2)m b 5 _ 16 7 2 U ' 1J (2, 1, 2)m 6 



^ \JZ y/Z 

+g 4°' 11 (3,l i 2)m fe 3 _ ig ff- 1] (2,l,2)m ft 3 _ g4 7 2 (1, 1, 2)m 6 5 



V 2 

+4 4°' 2l (3,2,2)m fc 4 _ 4 I 2 (3,2,l)m b 5 _ g 7 2 (2, 2, l)m b 5 



IZ y/Z y/Z 

i6 4°' 11 (2 i 2 ! 2)m b 5 | i6 /j°' 2l (3,l,2)m fc 3 g ' 2] (3, 2, 2)m b 3 



1 z z z 

|2i /^ 2] (2,2,2)m b 3 6 gjuV | 4 g^3 1 M>^! 

z z z 

64 7 2 °' 1] (1, 2, 2)m b 5 ^ f| ' 1] (2, 2, l)m b 3 ^ T^l, 1, 2)m b 5 

z z z 

| 2 7 (3, 1, l)m b 3 g A (2, 2, l)m b 5 7 (4, 1, l)m fe 3 

z z z 

|l6 4°' 2] (3,2,l)m b 3 | 2 J (2,2,2)m fe 5 | 7 (3, 2, 2)m fe 5 

2 z ' z 
21 / (l,3,l)m h 5 | I 2 (3,2,2)m b 6 | 16 7 (1, 2, 2)m 6 5 

z z z 

g ffl 1] (3,l,2)m 6 3 3 /l M (4,MK | . tT'&^K 5 

z z z 

4 J 1 (3, 2, l)m fc 5 2 7 2 (3, 1, 2)m b 5 ^ fj ' 1] (2, 2, l)m b 4 
z z z 

|g / 1 (2,l,2)m b 5 | ^2(3,2,1^ l J 2 °' 2] (3,2,2)m fe 4 

z z z 

|12 ffl 1] (3,l,2)m b 4 32 / 1 (2,l,l)m fc 4 | ^ I 2 °' 1] (2, 2, 2)m b 5 
z z z 

12 / 2 °' 11 (3,l ) 2)m b 4 | 32 / 2 (2,2,l)m b 5 I (3,l,2)m b 4 
z z z 

|20 / (2,2,l)m b 4 g lj°' 1] (3,2,l)m b 4 | g ffl 1] (2, 2, 2)m fe 4 

z z z 

g 7 2 (4, 1, l)m 6 4 /r'(3,2,2)m & 4 | ig if' 2 ' (2, 2, 2)m fe 4 
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, J 2 (2, 3, l)m fe 6 16 I (l, 2, l)m b 4 16 /j ' 1 ^, 2, 2)m fe 5 



+4. 



z z z 

7o(2, 1, 2)m b 4 ^ I 2 (l, 2, l)m b 5 ^ /f'^l, 2, 2)m fe 4 
z z z 

i6 4°' 11 (2,l i 2)m b 4 | 32 / 2 (l,2,2)m» 6 | 1 lj°' 2l (3, 2, 2)m b 4 

z z 



-12 



+ 4 

z z z 

I (l,4,l)m b 6 | 4 g^3 1 M>V | 3 / 2 (4,l,l)m fc 4 

z z z 3 / 2 

, /i°' 1] (2,2,2)m b 5 o / 1 (4,l,l)m b 4 h{Z,2,l)m h 5 

' 5 -3/2 + -3/2 



Z 3/2 " ^3/2 " ^3/2 

Tf^MjrV lj°' 1] (2,2,l)m b 4 / (2,2,2)m b 5 

£3/2 ^3/2 ^3/2 

A. 

^3/2 



-12 

•y^i^ 2°/^ ^o/^ 

ig /r'(3,2,2)m fc 3 , ff' 1] (l,4,l)m b 5 lj°' 2] (3, 2, l)m b 3 

+b ^72 + 48 ^72 + 16 ^ 

fi Ji(l,2,lK 5 ff ,1] (2,3,l)m 6 4 lj°' 2] (3,2,2)m fc 4 

z 3/2 8 ^3/2 4 ^3/2 

/{°' 11 (3,2,l)m fc 4 o/o / (4,l,l)m b 3 I 2 (3, 1, 2)m b 5 

-3/2 +<3/^ ^3/2 +^ ^3/2 



+4 



^3/2 
^72" 



/ 2 °' 2] (3,2,2)m b 4 A J^ 1] (2,2,2)m b 5 lj°' 1] (l, 2, 2)m b 5 

— 48 T^i b4 



-2 



2-3/2 £"V2 ^ 

/{°' 11 (l,2,2)m h 5 I (3,2,l)m b 4 I (2,2,l)m fe 

z 3/2 z 3/2 z 3/2 

I 2 (3,2,l)m fc 5 ig / (l,2,2)m fc 5 Qo ^(l, 2, 2)m b 6 

^3/2 + iD ^3/2 + ^ ^3/2 

J 1 (2,3 ! l)m fe 6 7 2 (l,2,2)m fc 6 l{°' 1] (4, 1, l)m fc 

^3/2 + ^3/2 ° z 3/2 



z 3/2 1 — z 3/2 ' z 3/2 

7 2 (l, 3,l)m b 6 7 1 (2,2,l)m b 5 I 2 °' 1] (3, 2, 2)m b 

.3/2 + ^3/2 + 4 ^3/2 

1 Zi , 

^3/2 1 ^3/2 

/x(3, 1, 2)m fe 5 I 2 °' 1] (2, 3, l)m b 5 I 2 (2, 2, 2)m b 6 

z 3/2 " ^3/2 8 z 3/2 

l0 , /o(M,lK , /i(3,2,2)m b 6 /io /j°' 2] (2,2,2)m fe 4 

+Z4 ^ + ^3/2 +4S z 3/2 

7 



+128^^ + 32-^ 

y (2,2,2k 4 o y(3,2,K , c< , Z (l,3,l)m» 5 
^72 8 ^72 + 64 



^3/2 

/ 2 (3,2,2)m b 6 ^(VMjW 1 nfi g' 1] (l,4,l)m & 6 

z 3/2 +4S z 3/2 +y ° z 3/2 



.3/2 



i0 /r'(3,2,2)m fe 4 io / 2 °' 1] (3,2,l)m b 4 J (2, 2, 2)m b 5 

+Z z 3/2 iZ z 3/2 +D .2 



-2 



£(3,2, l)m b 5 4g lj°' 1] (2,2,2)m b 5 ^ /j°' 1] (3,2, l)m b 4 



z A 
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| ^(3,1,2)77^ I 2 (3,2,2)m b 6 | d 7 (2, 3, l)m & 5 

z 2 z 2 z 2 

i0 / 2 (3,2,l)m b 5 /o / (3,2,2)m b 5 AjUW 
+ 2 p 1/2 ~ 2 + 48 

Jl?> 1] (3,2,2)m b * 

z 2 



9 li(4,l,lK ^ 72(1,4,1^ 

" + 3 i 2 " 48 ^ 

|l2 /{°' 2] (3,2,2)m b 4 6 / 2 °' 1] (3,2,2)m fc 5 | ^ 7 2 (2, 2, 2)m, 6 

z 2 z 2 z 2 

_ 24 7 (l,4,l)m & 6 _ g J 1 (2,2,2)m b 6 + , 

£ 2 Z 2 A 

|06 /{°' 1] (l,4,l)m b 6 | i /i°' 1] (3 i 2 ! 2)m b 5 ^ /{°' 11 (2, 3, l)m b 5 
z 2 z 2 
/!(3,2,2)m b 6 

-v2 



2 5 

/ 2 (2,3,l)m ft 6 

^ 2 



z 2 

6 



,2,2)m b 6 | 12£ /i(l,3,lK | ^ A (1,2,2)7/^ 
z 2 z 2 z 2 

i0 / 1 (3,2,l)m b 5 7i(l^l)W , 10 /i(2,2,2)m b 6 

+ 2 ^ 48 ^5/2 + 12 ^ 

I (3,2,2)m fe 5 Jj ' 11 (3 ) 2,2)m fc 5 I 2 (3 ; 2,2)m b 6 

+ / ^5/2 ^5/2 + ^5/2 

J 1 (3,2,2)m b 6 g J 1 (2,3,l)m & 6 J 1 (3,2,2)m fe 6 
+4 7 (2, 1, 2)m b 4 + 2 7 (3, 1, l)m b 3 + 7 (3, 2, l)m fe 4 

o 7 2 (2,l,2)m h 5 QO / 2 (2,l,l)m fc 4 7 (2, 1, 2)m fe 4 
— o — h oz — 4 — 

V 2 V 2 
ffl 1] (2,2,2)m b 4 J 2 °' 11 (3,l,2)m fe 4 7 2 (l,l,2)m ft 5 
/7 /7 /7 

1fl 71^(3, 1,2K 3 c ffl 1] (3,l,2)m 6 3 o /j°' 11 (2,l,2)m fc - 

b 7t b 7% v 8 7t 

V 2 v z 

i0 / 2 °' 1] (3,2,l)m b 4 ig / 2 °' 11 (2,2,2)m fe 5 7 (2, 2, l)m b 4 
+ 8 7r h ib 7r h 8 7r 

i1c 7 2 °- 1] (2,l,2)m 6 4 o / 2 °' 11 (3,2,2)m fe 5 7 2 (3, 2, l)m & 5 

V 2 V 2 V 2 

l / 2 °' 2l (3,2,2)m b 4 | g / 2 (2,2,l)m fe 5 ff' 1] (3, 2, 2)m fe 4 

v/i y/z y/z 

l / 2 °' 11 (3,2,l)m b 4 | l / 2 °' 2] (3 ! 2,2)m b 4 16 7 (1, 2, 2)m fe 5 

1c ffl 1] (2,l,2)m ft 4 
—lb 



z 

)m b 4 7 2 (3,2,2)- 6 



)m b 6 1c ffl 2] (3,l,2)m 6 3 
— ±u h lb 

z z z 

g 7l 0,11 (3,2,l)m & 4 | 3 / 2 (4 ; l ; l)m b 4 ^ 7 2 (2, 2, l)m fe 5 

z z z 
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4°' 1] (3,l,2)m fe 4 | i6 /f' 1] (2,2,l)m b 4 ^ h(3, 2, l)m b 5 
z z z 

I (2,2,2)m b 5 + J (3,l,2)m b 4 _ ^ lj°' 1] (2, 2, 2)m b 5 

z z z 

I 2 (3,2,l)m b 5 | ^f'ta.UK , 2 /r ] (3,2 ; 2)m fc 3 
z z z 

61 l,2)m fc 5 ^ Jj°' 11 (2,2,2)m fe 5 g 4°' 1] (3,2, l)m b 3 

z z z 

| 4 (3, 2, 2) W g (2, 2, V 16 If' 21 (3, 2, l)m» 
z z z 



+12 

-2 

-4 



+12 



z z z 

/ (l,4,l)m fe 6 o / (2,l,2)m 6 4 , 4°' 1] (3, 2, 2)m b 
o 1 

z z z 

I (4,l,l)m b 3 | g / 1 (2 ! l,2)m b 5 | g I 2 (2, 3, l)m fe 6 

z z z 



+3/2 

2 J (3, 2, 2)m fe 5 ig J 2 °' 2] (2, 2, 2)m b 4 g jfj (2, 2, 2)m b 4 

z z z 

|6 ffl 1] (4,l,l)m» 3 | ^ /f^l^m, 5 32 / 1 (2,l,l)m b 4 



2 Z Z 

I 2 (3, 1, 2)m b 5 ^ J (2, 2, l)m fc 4 ^ I 2 (l, 2, 2)m h < 

z z z 



+2 

Z 6 6 

|2 /j°' 1] (3,2,2)m b 5 | 2i / (l,3,l)m b 5 g A (2, 2, l)m b 5 

z z z 

|i2 (C , 1 )(3,l,2)m fe 4 | 6i / 2 (l,2,l)m b 5 | ig /|°' 11 (2, 2, 2)m fe 5 

z z z 3 / 2 

ioo / 1 (2,2,l)m 6 5 o / 1 (2,3,l)m 6 6 lj°' 2] (3, 2, l)m b 3 

+ 32 ^72 8 ^72 + 16 ^72 

Qfi |^(WH! on / 2 (l,2,2)m b 6 f^UjK 

" 96 ^72 32 ^72 64 ^ 

i , /j ' 1] (3 i 2,l)m b 4 ig / 2 °' 11 (2,3,l)m fc 5 Qo £(1, 2, 2)m b 6 

+4 ^3/2 + iD ^3/2 + ^ ^3/2 

I 2 (3,2,l)m b 5 _ J (l,3,l)m fe 5 I 2 (3,2,2)m b 6 

+ ^3/2 d ^3/2 + z 3/2 

h(3,2,l)m b 5 I 2 (3,l,2)m b 5 / 2 °' 1] (3, 2, 2)m fe 5 

+ z 3/2 z 3/2 z 3/2 

/}°' 21 (3, 2, 2)m b 4 4 ' 11 (3, 2, 2)m b 4 jf^ (3, 2, l)m fe 4 

z 3/2 z 3/2 + ^3/2 

/r ] (3,2,l)m b 3 Z 2 (4,l,l)m» 4 /j°' 11 (2 i 2, 2)m fe 5 

^3/2 d ^3/2 +iD ^3/2 

ig /{°' 1] (2,2,l)m b 4 o / (2,2,l)m fc 4 I 2 °' 2] (3, 2, 2)m b 4 

iD ^3/2 +S ^3/2 12 ^3/2 
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+g I 2 (2,2,2)m b 6 _ i2g I 2 (l,3,l)m fe 6 _ g A (4, 1, l)m b 4 



^3/2 



^3/2 



x3/2 



/J ' 11 (2, 3, l)rV _ 6 zf^&l.lK 3 _ 12 (3, 1, 2)m b 4 



^3/2 



^3/2 



^3/2 



+4 I (2,2,2)m fe 5 _ 2 J 1 (3,l,2)m b 5 + ^ /J ' 21 (2, 2, 2)m b 4 



^3/2 



^3/2 



^3/2 



^(l^,!^ 5 J (3,2,l)m fe 4 / (l,4,l)m fe 6 

^3/2 ^3/2 ^3/2 

_ 4g I 2 (l,4,l)m fc 7 + J 1 (3,2 ! 2)m fe 6 _ ^ (3, 2, l)m 6 4 



^3/2 

| 1 /{°' 1] (3 ! 2,2)m b 5 i6 /r i (2,3,l)m b 5 | d / (2, 3, l)m„ 5 



^3/2 



+3 J 1 (4 ! l,l)m b 4 + i2g Ji(l,3,l)m b 6 _ 2 I 2 (3,2,l)m fe 5 
z 2 z 2 z 2 

g / 1 (2,2,2)m b 6 | Qg ^kM^ljm^ g / 2 (2, 3, l)m 6 6 



|32 A(l,2,2)m 6 6 J 1 (3,2,2)m 6 6 | g A(3, 1, 2)m 6 5 

z 2 z 2 z 2 

|l2 /{°' 2l (3,2,2)m fe 4 | J (3,2,2)m b 5 g ^(3, 2, l)m b 5 

z 2 z 2 z 2 

^ / 2 °' 1] (3,2,2)m b 5 ^ /j°' 11 (2,2,2)m b 5 I 2 (3,2,2)m fc 6 

+0 48 1 



i2 / 2 (2,2 ; 2)m fe 6 | lg / 2 (l,4,l)m fe 7 | ^ I^l, 4, l)m b 7 

^,2 z 2 z 2 

J 2 (3, 2, 2)m 5 6 _ _ 4g 4, l)m 6 7 A(3, 2, 2)m b 6 



z b/2 - ^5/2 ^5/2 

_ g ff' 1] (3,2,2)m & 5 + g /i(2,3,l)m b 6 + ^ h (2, 2, 2)m fe 6 



+2 



^5/2 ^5/2 

/i(3,2,lK , h(3,2,2)m b 6 



y 5/2 



+ 



^5/2 

-2/ (3,l,l)m 6 3 



-Jo(3,2,l)m b 4 -4/ (2,l,2)m 6 4 



32 7q(1, 2, 2)m, 5 4 (3, 1, 2)m fc 3 ^ g ' 1] (3, 2, l)m fe 3 



/2 V 2 V 2 

+4 ff' 2l (3,2,2)m fe 3 + g J (2, 3, l)m 6 5 J (3, 2, 2)m fe 5 



z ^z yjz 

|g / (2,2,l)m b 4 32 /f' 1] (2,l,2)m b 3 | g I (4, 1, l)m fe 3 



z 



z 



g lj°' 2l (3,l,2)m 6 2 _ 2 / (3,l,2)m b 4 _ ^ 7^(2, 2, 2)m fe 4 



r[0,l] 
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where 



_ 16 / (2,l,l)m b 3 + g4 7 (l,l,2)m 6 4 + g ^(2, 1, 2)m 5 4 



'Z y/Z ^Z 

16 ff' 11 (2,3,l)m ft 4 g 1^(3,2, l)m b 2 | g 7 (3, 1, 2)m fc 4 

z z z 

| 6 Jp(3, 2, l)m b 4 g if > 1] (3, 1, 2)m b 3 ^ 7 (1, 4, l)m fe 6 

z z z 

^{3,2,2^ | 16 / (2,l,2)m & 4 | 16 7j°' 1] (3, 1, l)m b 2 

z z z 

| 1 /f |11 (3,2,l)m 6 3 | g / (3,l,l)m b 3 32 7 (2, 1, l)m 6 3 

z z z 

1Q2 7q(1, 2, l)m b 4 g 7 (2, 2, l)m fe 4 | g 7 (2, 2, 2)m b 5 



i6Q I (l,3,l)m b 5 + 2 / (3,l,2)m fe 4 _ ^ 7 (3, 2, l)m b 4 



^3/2 



^3/2 



+g J (2,3,l)m fc 5 + 7 (3,2,2)m & 5 _ g 7„(3, 1, l)m b 3 



+8 



^3/2 ^3/2 ^3/2 

4°' 1] (3,2,l)m & 3 o / (3,2,l)m b 4 



^3/2 

r[0,l]/o o o\... I 



+ 4 7 (3,2,l)m 6 4 



-2 7<T J (3, 2, 2)m b 4 + 8 7 (3, 1, l)m b 3 - 12 /f' 11 (3, 1, 2)m fe 3 
+4/ (2,2,2)m b 5 + 16/o(2,l,2)m fe 4 



r|^(a,6,c) 



(2m 6 )^ (jaVfja^ * 



d j 



2\3 



d(T?yd(Ti) 



Tl) MYi n (a,b,c) 
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Figure 3: The dependence of the form factor / on Borel parameters M\ and Mf for 
B c -> D 8 l+r jvv. 
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Figure 4: The same as Fig. 3, but for /_ . 
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fr 




Figure 5: The dependence of the form factor f T on Borel parameters M\ and Mf for 
B c ^D s l + l-. 
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Figure 6: The dependence of the form factor / on Borel parameters M\ and Mf for 
B c -> Dl + l~/vV. 
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Figure 7: The same as Fig. 6, but for /_ . 
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Figure 8: The dependence of the form factor f T on Borel parameters M\ and Mf for 
B c ^Dl + l-. 
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